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Abstract 

In this paper, we search for the regions of the phenomenological minimal supersymmetric 
standard model (pMSSM) parameter space where one can expect to have moderate Higgs mixing 
angle (a) with relatively light (up to 600 GeV) additional Higgses after satisfying the current 
LHC data. We perform a global fit analysis using most updated data (till December 2014) from 
the LHC and Tevatron experiments. The constraints coming from the precision measurements 
of the rare b-decays Bg —>■ and 6 —)• sy are also considered. We find that low (^ 350) 

and high tan/3 (> 25) regions are disfavored by the combined effect of the global analysis and 
flavour data. However, regions with Higgs mixing angle a ^ 0.1 - 0.8 are still allowed by the 
current data. We then study the existing direct search bounds on the heavy scalar/pseudoscalar 
(H/A) and charged Higgs boson (H^) masses and branchings at the LHC. It has been found 
that regions with low to moderate values of tan /3 with light additional Higgses (mass < 600 
GeV) are unconstrained by the data, while the regions with tan /3 > 20 are excluded considering 
the direct search bounds by the LHC-8 data. The possibility to probe the region with tan /3 < 
20 at the high luminosity run of LHC are also discussed, giving special attention to the H —hh, 
H/A —>■ tt and H/A —>■ r+T“ decay modes. 
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1 Introduction 


The ATLAS and CMS collaborations have now confirmed the discovery of a charge neutral spin-0 
particle with mass close to 125 GeV using the 7-1-8 TeV data at the LHC [1,2]. The CDF and 
D0 collaborations at the Tevatron have also reported the evidence of a SM-like Higgs boson with 
mass around 125 GeV [3] with Higgs boson decaying to bb. Measurement of it’s various couplings 
with the Standard Model (SM) particles so far seems to be consistent with the predictions of the 
SM Higgs boson. The primary goal of the run-H of the LHC with increased center-of-mass energy 
and integrated luminosity would be the more detailed study of various properties of the observed 
Higgs boson. The production modes of the observed Higgs boson that are analyzed at the LHC 
are the gluon-gluon fusion (ggF), vector boson fusion (VBF), associated production with a W or 
Z bosons, and associated production with a top-antitop pair, while the decay modes of the Higgs 
boson that are studied by the ATLAS and CMS collaborations^ are 77 [6,7], WW [8,9], ZZ [10,11], 
bb [12,13] and r+T“ [14,15]. Note that, the di-photon and ZZ modes have been used to precisely 
measure the mass of the observed Higgs boson [16,17]. Even though various coupling measurements 
at the LHC have already ruled out large deviations from the SM expectations, still the presence 
of uncertainties in coupling measurements do not rule out the possibility of having non-standard 
couplings of the observed Higgs boson. In fact, any small but statistically significant deviation from 
the SM expectations can be thought of as the first indication of new physics. 

Supersymmetry (SUSY) [18-20] has been one of the most popular extensions of the SM, however 
a SUSY signature is yet to be observed at the LHC. Non-observation of sparticles at the LHC have 
placed severe constraints on various superparticle masses and couplings [21,22]. The Higgs sector of 
the CP conserving minimal supersymmetric standard model (MSSM) contains two CP-even neutral 
Higgs bosons h and H, one CP-odd neutral Higgs boson A, and two charged Higgs bosons H^. One 
can identify the observed 125 GeV Higgs boson with the lightest MSSM Higgs boson h. At the 
tree level, the Higgs sector of MSSM is described by two parameters: pseudoscalar mass Ma and 
tan/I, where tan/3 is the ratio of the vacuum expectation values (vevs) of the two Higgs doublets 
Hu and The mixing angle a between the neutral components of the two Higgs doublets can 
be determined in terms of Ma and tan/3 at the tree level. However, radiative corrections to the 
Higgs boson mass matrix involving various SUSY parameters can modify the tree level value of a 
significantly. Moreover, the couplings of the lightest Higgs boson h with the SM particles depend 
on a and /3, and thus a global fit analysis considering various Higgs coupling measurements at the 
LHG and Tevatron experiments can, in principle, constrain the MSSM parameter space. 

Soon after the discovery of the 125 GeV Higgs boson, several analyses have been performed, for 
example, in the context of MSSM [23-38], general 2-Higgs doublet models [39-45], next-to minimal 
supersymmetric standard model (NMSSM) [46], effective theory framework [47-53]. However, in 

^Recently, search for the SM-like Higgs boson decaying to fJ.~ and e^e“ are also performed by both the ATLAS 
and CMS collaborations with 7-1-8 TeV LHC data [4,5]. 
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the last few months, both the ATLAS and CMS collaborations have updated some of their analyses 
on Higgs signal strength measurements. For example, ATLAS has updated their results for Higgs 
decaying to 77 [ 6 ], WW [ 8 ], ZZ [10], bb [12] channels, while CMS has published their new result 
for Higgs decaying to 77 [7]. In most of their analyses, the measure of uncertainty associated to 
various Higgs couplings have been now reduced by a sizable amount. A global analysis with the 
most updated (up to December 2014) results on various Higgs coupling measurements at the LHC 
and Tevatron would be extremely useful to probe the MSSM parameter space. 

In order to probe the Higgs sector of the MSSM, the discovery or exclusion of additional Higgs 
boson is extremely crucial at the LHC. In addition to several extensive studies of the sparticles 
of the MSSM, both the ATLAS and CMS collaborations have also performed dedicated searches 
of these additional Higgses in various possible hnal state signatures at the LHC. Searches for the 
heavy Higgses (both H and A) are performed at the LHC when they are produced via the gluon- 
gluon fusion and b-associated production processes with their decay to a pair of rs [54,55]. The 
ATLAS and CMS collaborations have also searched for the heavy resonances with the following 
decay modes: 77 —^ 77 [56], 77 —^ h/i —^ 6677 [57] and 77 —)■ /i/i —)• bbbb [58]. Heavy SM-like Higgs 
bosons decaying to pair of W bosons with final state consisting of one lepton, two jets and missing 
transverse energy has been studied by the CMS collaboration [59]. Heavy pseudoscalar Higgs boson 
(A) decaying to a Z boson and a light Higgs boson h has also been searched by both the ATLAS and 
CMS collaborations at 19.7 fb“^ luminosity at the 8 TeV run of LHC [60,61]. In addition, search 
for the charged Higgs bosons (H^) decaying to TUr, cs and tb are also performed at the LHC when 
H± are produced via tt or in association with top quarks [62,66]. So far, the ATLAS and CMS 
data have not revealed any clear signature of the signal and thus they put model-independent 95% 
C.L. upper limits on the production cross section times branching ratios for different production 
processes and decay modes. 

In the MSSM, the couplings of h with the SM electro-weak gauge bosons are proportional to 
sin(/3 — a). Now, the measured values of the couplings of the observed 125 GeV Higgs boson with 
the SM WjZ bosons are quite consistent with the SM expectations, which thus restricts the value 
of sin(/7 — a) near to unity. Besides, non observation of additional Higgses (77, A and 77^) at the 
LHC implies that their masses are possibly well above the electro-weak scale. One can satisfy both 
the above-mentioned observations in the MSSM by choosing a ~ 0 and /I ~ f and/or Ma » Mz, 
which is generally known as the decoupling limit of the MSSM [20]. In this paper, however, we 
would like to address the following question: Is there any MSSM parameter space still allowed 
by the current LHC data where one can expect to have moderate Higgs mixing angle (a) with 
relatively light (say, hundred to few hundred GeV) additional Higgses ? In other words, we are 
looking for a feasible MSSM parameter space where the non-negligible Higgs mixing still exists, 
and simultaneously we have light additional scalar particles which can be probed at the run-H of 
LHG. Our strategy can summarized as follows: 
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• We first perform a global analysis by scanning the relevant parameters of the MSSM 

Higgs sector and incorporating the updated Higgs signal strengths from the ATLAS, CMS 
and Tevatron experiments. The constraints coming from the precision measurements of rare 
decays like Bg —)■ and b ^ s'j are also considered. 

• We then impose the LHC direct search bounds on the heavy scalar/pseudoscalar Higgs (H/A) 
and charged Higgs boson (H^) masses and various branchings in the allowed parameter space. 

• Finally, we study the possibility to probe the remaining parameter space in the high luminosity 
run of LHC. 

In Sec. 2, we discuss the detailed prescription of our global fit analysis, followed by a brief 
outline of our parameter space scan. The favorable parameter space obtained after the global 
analysis is discussed in Sec. 2.2. In Sec. 3, we present the current limits on production cross section 
times branching ratios obtained from the direct search of both neutral and charged Higgs bosons 
at the LHC. Future limits for the high luminosity run of LHC are discussed in Sec. 4. Finally in 
Sec. 5 we summarize our results. 


2 Global analysis and available pMSSM parameter space 

In this section, we discuss the details of our global fit analysis. Current bounds on various Higgs 
coupling measurements by the ATLAS and CMS collaborations at the LHC, and also by the CDF 
and D0 collaborations at the Tevatron experiment are considered as the inputs to our global 
analysis. We also consider two important flavour physics constraints, namely Br (6 —)• sj) and 
Br(Hs —)■ ) [67]. Finally, we discuss the features of the available parameter space satisfying 

the updated Higgs and flavour physics data. 

2.1 Experimental inputs and Global fit analysis 

Both the ATLAS and CMS collaborations have published the results of the 125 GeV Higgs boson 
searches combining the 7+8 TeV data at the end of 8 TeV run of LHC [6-17]. Besides, the results 
of the SM Higgs boson search by the CDF and D0 collaborations at the Tevatron experiment are 
also available in the literature [3,68]. In our global analysis, we consider the most updated Higgs 
data obtained from the LHC and Tevatron experiments. The Higgs bosons are produced at the 
LHC mainly via the gluon-gluon fusion (ggF) process. However, there exist other sub-dominant 
production mechanisms e.g., vector boson fusion (VBF), associated prodnction with ^WjZ boson 
(Vh), associated production with a pair of top quarks [tth). The decay modes of the Higgs boson 
which are analyzed by the ATLAS and CMS collaborations are /i —^ 77 [ 6 , 7], /i —)■ WW* [8,9], 
h — 7 - ZZ* [10,11], h ^ bb [12,13], and h —?■ t^t~ [14,15], while the CDF and D0 collaborations 
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have analyzed the 77 , WW* and bb decay modes of the Higgs boson [3,68]. The experimental 
findings are usually presented in terms of the signal strength variable {p), which is defined as the 
ratio of the production cross section (a) times the branching ratio (Br) to a specific decay mode for 
a given new physics model normalized to the SM prediction. For example, when the Higgs boson 
is produced via gluon-gluon fusion process and it decays to a generic final state XX ( 77 , WW*, 
ZZ*, bb and then one can define the signal strength variable p, assuming narrow-width 

approximation, as: 


PggPiXX) 


T{h gg) Br{h XX) 

r(/isM gg) ^ Br{hsM xx) ’ 


( 1 ) 


where /i is a observed 125 GeV Higgs boson and hgM is the SM Higgs boson. Similarly, if the Higgs 
boson is produced via VBF fusion process and it decays to XX, then one can define. 


g-VBF/VHi^^) 


r{h WW) Br{h XX) 
r(/i5M WW) Br{hsM XX) 


( 2 ) 


From Eq. 1 and 2 it is evident that the signal strengths are functions of the partial decay widths 
and the total decay width of the Higgs boson. The presence of new particles/interactions in the 
new physics models leads to modifications in the partial/total decay widths, and thereby changes in 
the value of signal strength variable. Thus, precise measurements of these signal strength variables 
are extremely crucial as a small but statistically significant deviation from the SM expectation will 
hint at possible signatures of the new physics. Note that, in the timeline of the Moriond 2014 [69] 
workshop to the end of December 2014, significant changes have been observed in the measurement 
of various Higgs signal strengths. For example, ATLAS data for di-photon signal strength has 
changed from 1 . 57 ]’jg' 2 g [69] to 1.17 ± 0.27 [ 6 ], while the same from the CMS has changed from 
0.77 ±0.27 [69] to l-UtoH [7]. Besides, both statistical and systematics uncertainties associated 
to some of these signal strengths have been reduced after the combination of (7±8) TeV LHC data. 
Near the end of last year, the CDF and D0 collaborations have also updated the constraints of 
125 GeV Higgs boson couplings to fermions and vector bosons. In this work, we perform a global 
analysis considering all these most updated Higgs signal strengths as of the end of 2014. The 
available signal strength variables^, for different production and decay modes of the Higgs boson 
at the LHC are summarized in Table 1 to Table 5 while Tevatron Higgs data is shown in Table 6 . 

Let us now discuss the details of the MSSM parameter space scan. We do not consider any 
specific SUSY breaking scenario, rather we focus on the generic phenomenological MSSM (pMSSM) 
model with 19 free parameters and perform a random scan for approximately 100 million points. 
Those parameters that are relevant to the MSSM Higgs sector, namely pseudo-scalar mass pa¬ 
rameter Ma, the ratio of the vacuum expectation values of two Higgs doublets tan/3, higgsino 

^For global fit analysis, we consider the signal strengths for different production modes as presented for individual 
decay channels. In other words, we do not consider the inclusive results for a given decay mode. 
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Channel 

Signal strength (/i) 

Production mode 

ATLAS CMS 

ggF VBF Vh 

Kagh) 

H{VBF) 

n{Wh) 

1.32 + 0.38 

0.8 ±0.7 1.58l[3;^^ 

1.0 ±1.6 -O.ielJ;^^ 

0 1+3-7 
^•-‘--0.1 

100 % 

100 % 

100 % 

100 % 


Table 1: Signal strengths of /i —)■ 77 channel as recorded by the ATLAS [ 6 ] and CMS [7] collaborations 
after 7+8 TeV run of LHC with 25 fb“^ of luminosity. The amount of contribution to a given channel from 
each production modes are shown in Column 4-6. The total for the 77 channel with respect to the SM 
is = 0.85 (ATLAS) + 1.868 (CMS) = 2.718. 


Channel 

Signal strength (/i) 

Production mode 


ATLAS 

CMS 

ggF 

VBF 

Vh 

n{ggh ± bbh ± tth) 
g{VBF + Vh) 

1 66+3-3^ 
r-OO_0.44 

0 26+^-®^ 
^•^0-0.94 

0.80^°;^® 

1 '7+2.2 
-2.1 

100 % 

60% 

40% 


Table 2; Signal strengths of /i —>■ ZZ* channel as recorded by the ATLAS [10] and CMS [11] collaborations 
after 7+8 TeV run of LHC with 25 fb“^ of luminosity. The amount of contribution to a given channel from 
each production modes are shown in Column 4-6. The total for the ZZ* channel with respect to the SM 
is = 2.493 (ATLAS) + 0.3 (CMS) = 2.793. 


Channel 

Signal strength (g) 

Production mode 


ATLAS 

CMS 

ggF 

VBF 

Vh 

TiggF) 

1 02+3-^® 

- 

100 % 

- 

- 

g{VBF) 

1 07 + 0.53 
-*-•^'-0.45 

- 

- 

100 % 

- 

/r(0/l jet) 

- 

0 74+0-22 
"o- '"±- 0.20 

97% 

3% 

- 

g (VBF tag) 

- 

o-eolliie" 

17% 

83% 

- 

g (Vh tag {2l2u2j)) 

- 

0.39l}|7 

- 

- 

100 % 

g (Wh tag(3/3i/)) 

- 

0-561^:95 

- 

- 

100 % 


Table 3: Signal strengths of h ^ WW* channel as recorded by the ATLAS [ 8 ] and CMS [9] collaborations 
after 7+8 TeV run of LHC with 25 fb“^ of luminosity. The amount of contribution to a given channel from 
each production modes are shown in Column 4-6. The total y^ for the WW* channel with respect to the 
SM is = 0.366 (ATLAS) + 2.104 (CMS) = 2.470. 
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Channel 

Signal strength (//) 

Production mode 

ATLAS CMS 

ggF VBF Vh 

/u(Vh tag) 

1.0 ±0.5 

100% 


Table 4: Signal strengths oi h ^ hh channel as recorded by the ATLAS [12] and CMS [13] collaborations 
after 7+8 TeV run of LHC with 25 fb“^ of luminosity. The amount of contribution to a given channel from 
each production modes are shown in Column 4-6. The total for the bb channel with respect to the SM is 
= 1.50 (ATLAS) + 0.0 (CMS) = 1.5. 


Channel 

Signal strength (/r) 

Production mode 

ATLAS 

CMS 

ggF 

VBF 

Vh 

K 99 F) 

fi{VBF + Vh) 

1 qQ+l-45 

-‘-•^^-0.54 

- 

100 % 

60% 

40% 

9 (0-jet) 

- 

0.34 + 1.09 

96.9% 

1.0% 

2.1 

/i (1-jet) 

- 

1.07 + 0.46 

75.7% 

14% 

10.3 

/i (VBF tag) 

- 

0.94 ±0.41 

19.6 

80.4 

- 

/i (Vh tag) 

- 

-0.33 + 1.02 

- 

- 

100% 


Table 5: Signal strengths of —>■ r+r“ channel as recorded by the ATLAS [14] and CMS [15] collaborations 
after 7+8 TeV run of LHC with 25 fb“^ of luminosity. The amount of contribution to a given channel from 
each production modes are shown in Column 4-6. The total fbe channel with respect to the 

SM is = 0.857 (ATLAS) + 2.11 (CMS) = 2.967. 


Channel 

Signal strength (fi) 

Production mode 

Tevatron 

ggF VBF Vh 

77 ) 

^WW*) 

—)• bb) 

e.utlil 

1 t;n+0.69 
i.oy_0 72 

78% 5% 17% 

78% 5% 17% 

100% 


Table 6 : Signal strengths of ft. —>■ 77 , WW*, and bb channel as recorded by the CDF and 110 collaborations 
at the Tevatron with 10 fb“^ of luminosity at yT = 1.96 TeV [3,68]. The amount of contribution to a given 
channel from each production modes are shown in Column 4-6. The total for the above three modes with 
respect to the SM is = 3.296. 


mass parameter /i, the third generation sqnark trilinear couplings At and Ab (trilinear couplings 
of sleptons and first two generations squarks are set to zero), third generation squark soft mass 





















parameters Mq^, Mjj^ and are scanned in the following ranges: 


1 < tan /3 < 50, 100 GeV < Ma < 600 GeV, 

-8000 GeV < At, < 8000 GeV, 100 GeV < // < 8000 GeV, 

100 GeV < Mq 3 , Mu 3 < 8000 GeV, 100 GeV < Md 3 < 8000 GeV, (3) 

while we fix the following parameters since they have little impact on our analysis. 

Ml = 100 GeV, Ms = 2000 GeV, M 3 = 3000 GeV, 

-^Li 2 3 “ -^Ei 23 “ 3000 GeV, Mqj^ ^ — 3000 GeV, Muj^ ^ — -^Di 2 “ 3000 GeV, (4) 


where Mi^ 2,3 are the gaugino mass parameters, Ml; and Me; (i = 1,2,3) are the left and right 
handed slepton soft SUSY breaking mass parameters, and Mq., Mtj., Md. (i = 1,2) are the first 
two generation squark soft SUSY breaking mass parameters. 


We scan the third generation trilinear couplings {At and Af,) and soft masses (Mq 3 , Mu 3 , Md 3 ) 
over a wide range in order to obtain the lightest MSSM Higgs boson mass in the range of 125 ± 
3 GeV assuming 3 GeV uncertainty in Higgs mass calculation [70]. Since we are interested in the 
possibility of having light additional Higgses, we restrict Ma up to 600 GeV. From the above choices 
of the model parameters, it is evident that we do not consider the possibility of the decay of h to 
MSSM particles. We use SUSPECT (version 2.43) [71] to scan the MSSM parameter space and 
Superlso (version 3.4) [72] to calculate the flavour physics observables, while the branching ratios 
of the lightest Higgs boson are evaluated using HDECAY (version 6.41) [73]. 


Now, to combine the available information on different signal strength variables from the AT¬ 
LAS, CMS and Tevatron, and to compare with the MSSM expectations, we compute for all the 
scanned parameter space points, defined as below: 


2 _ ~ 

^ “A 


(5) 


where [it is the experimentally observed signal strength for a particular production/decay mode i, 
and Jit is the value predicted for the same channel for a chosen MSSM parameter space point with 
A/ij being the measure of the experimental error associated to that channel. The sum over i takes 
into account of all the experimentally measured production and decay modes of the Higgs boson. 


It is to be noted that, different production processes can, in principle, contribute to a particular 
experimental search channel, thereby while calculating the signal strengths for a parameter space 
point, contributions coming from different production processes need to be considered. Hollowing 
the procedure of Ref. [29], we implement this modihcation in our analysis as follows: 

= ( 6 ) 

where denotes the amount of contribution that can originate from the production mode j to 
the category/channel i with /ij being the signal strength corresponding to the MSSM parameter 
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space point. For example, the category /i(VBFtag), as introduced in Table 3, receives 17% and 
83% contributions from the ggF and VBF processes respectively [74], So, we calculate Jj-ggP and 
JlvBF, and then scale them with 0.17 and 0.83 (the T-s here) respectively to obtain the proper 
signal strength (jl) corresponding to VBFtag category. To obtain the contributions coming from 
different production processes to a given decay mode/category, we use LHC Higgs cross section 
working group report [74]. 

We consider altogether 28 data points (i.e., experimental inputs, see Table 1-6) combining the 
CMS and ATLAS and Tevatron Higgs data. We calculate foi' all the scanned parameter space 
points and find the minimum of x^- We call this minimum as the approximate minima (Xapprox)- 
In order to obtain the true x^ minimum (Xmin)’ vary the parameters around their approximated 
values i.e., values corresponding to the approximate x^ minimum. We present the parameter space 
that is available after the global analysis considering the la and 2cj intervals with x^ = Xmin + 2.3 
and x^ = Xmin + respectively in Ma — tan/3 plane [75]. The “best fit” value corresponds to 
Ma ~ 584 GeV and tan/3 ~ 36. To determine how a set of experimental data is well represented by 
any given model, one usually calculates the chi-square per degrees of freedom (d.o.f) i.e., x^/d.o.f. 
The minimum value of x^ obtained from the analysis for SM is 15.744 with x^/d.o.f = x^/28 = 
0.562, while for MSSM we obtain Xmin= 15.013 with x^/d.o.f = x^/20 = 0.75. 

2.2 Available MSSM parameter space 

Before we proceed to discuss our findings, let us review our methodology once more. We first 
perform a x^ analysis using a random scan of the parameters relevant to the MSSM Higgs sector, 
then we find the true Xmin get the “best-fit” values of those parameters. The points with x^ within 
2a of the true Xmin only considered for further analysis. We impose the two most stringent rare 
b-decay constraints, namely Br(b —)■ sy) and Br(Bs —)■ and allow 2a deviation^ [67], 

2.82 X 10“^ < Br(Bs ^ X^y) < 4.04 x 10“^ 

1.57 X 10"^ < Br(Bs ^ < 4.63 x 10"^ (7) 

In Fig.l, we show the parameter space in the Ma — tan/3 plane obtained from the global fit 
analysis and also satisfying the flavour physics constraints. Magenta (black) coloured triangle 
(circle) shaped points represent 2a (lu) allowed parameter space. The region with Ma < 350 GeV 
and tan/3 > 25 are excluded by the stringent Br(Bs —)• constraint, which is expected to 

dominate for regions with large tan/3 and relatively smaller Ma. However, for significantly larger 
Ma i.e., Ma > 425 GeV, the effect of this constraint is negligible. Besides, most of the points in 
the region with Ma < 350 GeV with tan /3 < 8 are excluded by the Br(b —)• sy) constraint. 

^The current measurements of these two b-observables are Br(i3s —>■ Xsj) = 3.43±0.22±0.21(theo.) and BrfBs —t 
U~) ~ 3.1±0.7±0.31(theo.) [67]. We follow Ref. [76] for the conservative estimates of the theoretical uncertainties 
associated to these two flavor observables. 
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Ma [GeV] 


Figure 1: Parameter space allowed in Ma - tan/d plane from our global fit analysis and also 
satisfying the flavour physics constraints on Br(b s^) and BrfBg —)• ). Magenta (black) 

coloured triangle (circle) shaped points represent 2a (la) allowed parameter space from global fits of 
125 GeV Higgs data after Run-I of LHC. In rest of our analysis, while presenting the direct search 
constraints and the future limits on the heavy Higgs masses and BRs, we would consider these 2a 
points which satisfies our global fit analysis and the updated flavour data. 



100 200 300 400 500 600 100 200 300 400 500 600 100 200 300 400 500 600 

Ma [GeVl Ma [GeV] Ma [GeV] 


Figure 2: Scatter plots in the (a) Ma — tan/3 plane without any flavor constraint (left panel), (b) 
Ma — tan/3 plane only after imposing Br(b sj) constraint (middle panel) and (c) Ma — tan/3 
plane after imposing only Br(Bs —)• constraint (right panel). In Fig. 1, we show the same 

correlation after the imposition of both Br(b s^) and Br(Bs —)• ) constraints. 
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Figure 3: Left panel (a) shows the allowed parameter space in tan/3 — a plane after the global fit 
analysis and also satisfying the flavour physics constraints. Colour conventions are same as Fig. 1. 
In the right panel (h), we display the same allowed parameter space but now is presented in the 
Ma — a plane. 


In order to display the interplay of two flavour physics constraints, in Fig. 2(a) we first show 
the parameter space obtained after the global analysis with Higgs mass constraint (122 GeV < 
M/j < 128 GeV) only, without imposing any flavor physics constraint. In the middle panel (6), we 
then show the same distribution but only after imposing the Br(b —)• sy) constraint. In right-most 
figure (Fig. 2c) we plot the same correlation with imposing Br(Bs —> pfipi~) constraint only. The 
effect of the Br(Bg —)• constraint in the low Ma and large tan/3 region, and the impact of 

Br(b —)• sy) in low Ma and low tan /3 is now clearly visible from these plots. Once we impose both 
the Higgs mass and flavour physics constraints, the available parameter space has already been 
shown in Fig. 1. Note that, in rest of our analysis, we name these 2(t allowed points collectively as 
the “scanned data set” and present the direct search constraints and the future limits on the heavy 
Higgs masses and couplings using this data set. 

In Fig. 3, we show the scatter plots in the (a) a — tan /3 and (6) a — Ma planes, where a is the 
Higgs mixing angle. But, before we proceed, let’s first briefly discuss the idea of the “alignment 
limit” (for details see Refs. [77-80]). In the MSSM, the couplings (at the tree level) of the GP-even 
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Higgs bosons {h,H) to SM gauge bosons are [20], 

9hvv = sin(/3 - a) gy 

gHVV = cos(/3 - a) gv, (8) 

while, couplings to the SM fermions are [20], 

9hdd = — sin a/ cos f3gf = (sin(/3 — a) — tan/?cos(/3 — a)) gj 
9huu = — cos a/ sm f3gf = (sin(/3 — a) + cot (3 cos(/3 — a)) gj 
9Hdd = — cos a/ cos f3gf = (cos(/3 — a) + tan j3 sin(/3 — a)) gf 
9Huu = - sin a/ sin Pgj = (cos(/3 - a) - cot jS sin(/3 - a)) gf, (9) 

where gy and gf are the corresponding SM couplings with gy = 2iMy/v and 5 / = iMf/v for a 
generic gauge boson V (V = W, Z) and a fermion / with v = 246 GeV. Now, the alignment limit is 
the limiting case when the lightest CP-even Higgs boson mimic the properties of the SM Higgs and 
the SM gauge bosons couple to the light CP-even Higgs boson only ( i.e., ghw ~ 1 ^-^icl gnvv ~ 0)- 
This limit can be easily achieved with the variation of the two quantities a and /3. For example, 
let’s assume a special case: {13 —a) ~ 7 r /2 for some specific choices of a and j3, then the couplings of 
the lightest CP even Higgs boson to the SM gauge bosons are SM-like, while heavier CP even Higgs 
boson couplings become highly suppressed. Note that, even though the heavier CP even Higgs has 
highly suppressed coupling to the SM gauge bosons, however it can have non-zero couplings to the 
SM fermions depending on the choice of tan/3 (see Eq. 9). Now, let’s discuss how we achieve the 
alignment limit from our parameter space scan. To do so, let’s first divide the region of interest 
into three parts, namely tan/3 < 5, 5 < tan/3 < 20, and tan/3 > 20. The region with tan/3 < 5 
corresponds to a > - 0.2 radian. Now, tan/3 = 5 implies /3 = 1.373 radian, which means (/3 — a) = 
1.573 radian or, around 90 degree with a = - 0.2 radian, so here we are near the alignment limit 
{P — a = vr/2) [77]. The region with tan,3 > 20 corresponds to very small a (< - 0.05 radian), and 
thereby {(3 — a) ~ tt/2 i.e., we again achieve the alignment limit. The intermediate regime with 
5 < tan /3 < 20 also satisfies the criteria of alignment. However, we see that in this alignment limit, 
Ma can be as light as 300 - 400 GeV with relatively large a (see Fig. 3(6)) satisfying current data. 
This scenario can be thought of as the “alignment without decoupling” scenario, as discussed in 
Ref. [78-80]. In fact, in order to make a more quantitative statement, In Fig. 4 we present the 
allowed parameter space in the {(3 — a) — tan/3 plane with different choices of Ma- We divide the 
entire Ma region into four parts: 200 < Ma < 300 GeV (red/square), 300 < Ma < 400 GeV 
(black/circle), 400 < Ma < 500 GeV (green/triangle), and 500 < Ma < 600 GeV (blue/cross), 
while the red horizontal line at {/3 — a) = 1.571 indicates the exact value at the alignment limit. 
From Fig. 3 and Fig. 4 it is clear that regions with light Ma (< 400 GeV) satisfying the alignment 
limit is perfectly allowed by the current data, one is thus not always forced to be in the decoupling 
limit to comply with LHC data. 
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Figure 4: The allowed parameter space in the (/3 — a) — tan/3 plane is presented with different 
choices of Ma- The entire Ma region is divided in four parts, namely 200 < Ma < 300 GeV 
(red/square), 300 < Ma < 400 GeV (black/circle), 400 < Ma < 500 GeV (green/triangle), and 
500 < Ma < 600 GeV (blue/cross). The red horizontal line at {fd — a) = 1.571 indicates the value 
of {a — P) at the alignment limit. We see that regions with light Ma (< 400 GeV) satisfying the 
alignment limit are perfectly allowed by the current data. 
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Figure 5: Signal strength correlations in (a) yzz — g-'y-f, (bjgzz — gr+T- and (c) gzz — gsb Magenta 

(black) coloured triangle (circle) shaped points represent 2a (la) allowed points of our “scanned data set”. 


Before we end this section, we would like to discuss the correlations of various Higgs signal 
strength variables. In Fig. 5, we show the scatter plots in gzz — h'j'y (left), gzz — hr+r- (middle) 
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and fizz — fJ-bb (right) planes. We find that the partial widths T{h —)• ZZ) and T{h —)• WW) remain 
almost unaltered with respect to the SM values, however r(/i —)• bb) and total decay width of Higgs 
{Ttot) have increased while r(/i —)■ gg) has decreased for most of the scanned data points. Now, we 
know that when the total decay width increases, the branching ratios of various sub-leading decay 
modes like 77 , gg, ZZ, WW decreases. Hence, the suppression in the fizz/fi-y-y can be thought 
of as the interplay of two effects: the increase in total decay width (or, decrease in /i —)• ZZ or 
/i —)• 77 branching ratios) and decrease in T{h —)• gg) with respect to the SM. Now, in order to 
understand the bb and rr correlations in Fig. 5(b) and Fig. 5(c), we need to discuss the couplings 
of the lightest CP-even Higgs boson with the b quark and r lepton. At the tree level, the Yukawa 
couplings of b and r are proportional to (see Eq. 9). However, loop corrections (in powers 
of as tan /3) involving heavier supersymmetric particles can significantly modify the b quark mass 
and it’s Yukawa coupling from its tree level predictions. In an effective Lagrangian approach, these 
effects are parametrized by the quantity A;, and one can write [81-84], 


T - - ( 1 

vsM \ 1 + Afo 


sma 
cos (3 


1 - 


Afe 


tan (3 tan a J 


\hh 


where vgM = {V^Gp) The bottom quark mass can now be written as, 

VbVi 


nib 


^/2 


(1 -|- Ab), 


( 10 ) 


( 11 ) 


where vi = vsm cos j3 and yb the bottom Yukawa coupling. In order to understand various terms 
of Eq. 10, in Eig. 6 (a), we show the distribution of the quantity. 
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1 + Af 



Af, 


tan /3 tan a 


( 12 ) 


with respect to A^ for all the 2a (magenta triangle) and Icr (black circles) allowed points. It is 
evident from the figure that for significant number of points. A;, is mostly positive and varies within 
10-15%. However, the effect of this variation of A^ on e is small. Thus, we understand that A;, is 
not playing a very significant role here, and so we then proceed to estimate the individual quantities 
that are involved in firr and fibb- Erom Eq. 1 and Eq. 2, we find that the quantities of interest are 
T{h —)• gg), T{h —>■ WW), BR{h —)• ZZ), BR{h —)• rr) and BR{h —)■ bb). In Eig. 6 (b), we display 
the variation of the ratio of partial widths of Higgs to bb and rr with respect to the SM values with 
the tree level coupling Eor most of the 2 it allowed points, we find 25 - 30% enhancement in bb 
and rr partial widths (see Eig. 6 (b)), however we check that F(/i —)• ZZ) remains unaltered. Thus, 
increase in F(/i —)• bb) and thereby increase in Ftot results in 20-25% modifications in BR{h —)• rr), 
however, being the dominant decay mode, change in BR{h —)• bb) is small (see Eig. 6 (c)). Thus, 
one can conclude that the larger spread in fi-rr as seen in Eig. 5(b) compared to Eig. 5(c) is coming 
from the interplay of the total Higgs decay width and individual Higgs branching ratios along with 
a mild dependence of A;,. 
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Figure 6: Left panel shows the distribution of the quantity e = x tan an a J with the variation 

of Afe. The variation of the ratios of the partial decay widths and branching ratios of the Higgs boson to 
the SM values for the bb and T+r“ modes with the tree level Higgs Yukawa couplings are shown in the 
middle and right panel respectively. For more details, see the text. 


3 Bounds on MSSM heavy Higgses from direct searches 

In the previous section, we first describe the global fit analysis, and then show that the signal 
strength measurements of the SM-like Higgs boson with different possible hnal state signatures 
do not exclude the possibility of having additional light Higgses (say < 600 GeV) with moderate 
values of Higgs mixing angle a. One can now ask whether these light MSSM Higgses are still 
allowed satisfying the direct search bounds at the LHC-8, which is precisely the goal of this section. 
Here we will impose the bounds set by the ATLAS and CMS collaborations on the masses and 
branching ratios of the neutral and charged Higgs bosons at the end of 8 TeV run of LHC. We 
calculate the production cross section of the neutral heavy Higgses (H and A) using SuShi (version 
1.4.1) [85] and charged Higgs using PYTHIA (version 6.4.28) [86]. The branching ratios of 

both the charged and neutral Higgses are evaluated using HDECAY (version 6.41) [73]. 

3.1 Neutral Higgs boson searches 

3.1.1 Search for H with 77 final states 

The di-photon invariant mass distribution plays an important role to discover the 125 GeV Higgs 
boson at the LHC. However, the sensitivity of this channel falls rapidly with the increase of the SM 
Higgs boson mass and become vanishingly small beyond 150 GeV [56]. In models with additional 
Higgses, di-photon mode can be a useful probe to search for these heavy resonances. The CMS 
collaboration has searched for the CP-even heavy MSSM Higgs H using the di-photon invariant 
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mass distribution with 19.7 fb“^ of data at the 8 TeV run of LHC. They assume that the Higgs 
is produced via gluon-gluon fusion process. Both narrow and wide width heavy resonances are 
investigated with widths ranging from 0.1 to several GeV and Higgs masses varying in the range of 
150 GeV to 850 GeV. No excess over the SM background has been found, and thus 95% C.L. upper 
bounds have been set on the production cross section times branching ratio in the above-mentioned 
Higgs mass range. In Fig. 7(a), we show the distribution of the quantity a x Br(H —)■ 77 ) for all 
the points of our scanned data set (i.e., 2a allowed points obtained after the global analysis with 
various Higgs signal strengths and flavour data, see Sec. 2). We then superimpose the bounds set 
by the CMS collaboration for two different choices of the Higgs decay width. The red solid line in 
Fig. 7(a) represents the case where Higgs decay width is 10% of the Higgs boson mass, while the 
blue dashed line represents the same but with fixed value of the Higgs decay width F = 0.1 GeV. 
From the figure, one can see that all the points of our scanned data set satisfy the CMS bounds. A 
further investigation reveals that the Higgs to 77 branching ratio for all the points corresponding 
to the scanned dataset varies between 10 “® — 10 “'^ for all the points, and thus makes the quantity 
a X Br(H —)■ 77 ) small enough to evade the CMS bound. 

3.1.2 Search for H with WW final states 

In the SM, the decay of Higgs boson to electro-weak gauge bosons (W/Z) provides interesting 
signatures at the LHC. In fact, the ZZ mode has a very good sensitivity to precisely measure 
the Higgs boson mass and its spin and parity [16,17]. The CMS collaboration has searched for a 
SM-like Higgs boson decaying to a pair of W bosons with (5 -|- 19.3) fb“^ data collected at ^/s = 
(7 -|- 8 ) TeV at the LHC [59]. Analyzing the data with lepton, jets plus missing transverse energy 
final state signature, CMS collaboration excludes SM-like Higgs bosons in the mass ranges 170-180 
GeV and 230-545 GeV at 95% G.L.. One can translate the bounds on the cross sections obtained 
from search of SM-like Higgs boson via this channel to the same corresponding to a model with 
additional Higgses. We calculate the product of production cross section and WW branching ratio, 
a X Br(H —)• WW), for our scanned points and then compare our findings with the CMS data. Note 
that we do not combine the 7 and 8 TeV data of LHC, rather we consider the exclusion limits of 
the 8 TeV data only. In Fig. 7(b), we superimpose the CMS exclusion limit on the 2a allowed 
points obtained from global fit analysis. We find that updated CMS bound on production cross 
section with Higgs decaying to pair of VFs is not sensitive enough to exclude the scanned parameter 
space. As we have already discussed, for almost all the points we are near the alignment limit (i.e. 
{13 — a) ~ |), which results into highly suppressed branching ratio of Higgs to the WW mode, 
Br(H —)• WW) ~ 10“^ —10“^. This is the reason why the quantity cr x Br(F7 —)• WW) is very small, 
and thus become insensitive to the CMS bound. We also find that even if we consider the (7-1-8) 
TeV combined CMS exclusion limit the parameter space will still remain allowed by the LHC data. 
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Figure 7: (a) Left: Scatter plot in Mh — [ax Br(H — )■ 77)] plane assuming gluon-gluon fusion 

production process. The solid red (blue dashed) line represents the observed upper limits on a x 
Br(H — )■ 77) at 95% C.L. by the CMS collaboration [56] using LHC-8 data with Th = 0.1 x Mh (0.1 
GeV). Magenta coloured triangle shaped points represent 2a allowed parameter space from global 
fits satisfying flavour constraints (see Fig 1). (b) Right: Scatter plot in Mh — [ct x Br(H —)• WW)] 
plane when the H is produced via ggF. The red solid line indicates the observed upper limits on 
a X Br(H —)• WW) at the 95% C.L. by the CMS collaboration [59] using LHC-8 data. 


3.1.3 Search for H with hh (bbbb and 6677 ) final states 

The Higgs pair production cross section in the SM is very small, around 10 fb at y/s = 8 TeV [57]. 
However, well-motivated BSM physics models predict the decay of narrow-width heavy resonances 
to a pair of 125 GeV Higgses, and thus one can expect enhancement in the 125 GeV Higgs pair 
production cross section. Search for these heavy resonances decaying to a pair of h, i.e. pp —?• V —?• 
hh with X being the new heavy resonance, is performed by the CMS collaboration in the mass range 
of 260 - 1100 GeV with 19.7 fb“^ of data at \/s = 8 TeV [57,58] The final state signatures that 
have been investigated by the CMS collaboration are: (z) bbbb where both the Higgses decaying to 
bh, and (ii) 6677 where one of the Higgs decays to bb while other decays to a pair of photons. The 
reconstruction of heavy resonance is possible in both the above-mentioned signatures, however 46 
final state is experimentally challenging while 6677 channel has less background contamination and 
very good di-photon mass resolution. The observations by the CMS collaboration are consistent 
with the SM, and so 95% exclusion limits on the production cross sections are placed for the 
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entire mass range of the heavy resonance. At the LHC, the MSSM heavy CP even Higgs (H) can 
be generously produced from the ggF process, which then decays to a pair of 125 GeV Higgses. 
So, the bounds on the production cross section of new heavy resonances can be directly used to 
probe/exclude certain mass range of the heavy Higgs boson. To do so, we calculate the quantities 
a{pp —)• H) X Br(H —?• hh —?• bbbb) and cr{pp — )• H) x Br(H —?• hh —)• bbyy) for all our scanned data 
set and in Fig. 8(a) and Fig. 8(b), we show the distributions of those two quantities respectively. 
The 95% C.L. upper limits (solid red line) imposed by the CMS collaboration on the respective 
distributions are also superimposed. 
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Figure 8: (a) Left: Scatter plot in Mh — [ax Br(H —)• hh —)• bbbb)] plane assuming gluon- 

gluon fusion production process for 2a allowed parameter space (represented by magenta points). 
The solid red line represents the observed upper limits on a x Br(H —)• hh —?• bbbb) at 95% 
C.L. by the CMS collaboration [58] using 8 TeV data, (b) Right: Scatter plot in Mh — [ax 
Br(H —)■ hh —)• bbyy)] for ggF production mode. The red solid line indicates the observed upper 
limits on a X Br(H —)> hh —)■ bbyy) at the 95% C.L. by the CMS collaboration [57] using LHC 8 
TeV data. 


The branching fraction of the MSSM heavy Higgs H to hh becomes sizable only for small tan /3 
(< 5) and low Ma [20]. Once Ma becomes > 350 GeV, the ti decay mode opens up and dominates 
in the rest of the parameter space. We find that so far CMS data can not exclude our 2cr allowed 
parameter space, however the regions close to CMS exclusion lines corresponds to smaller values of 
tan/I, and so with one/two orders of improved measurement of di-Higgs production cross section 
at the Run-H of LHC, one can probe/exclude such regions of the parameter space. The ATLAS 


19 


















collaboration has also searched for the heavy resonances (Kaluza-Klein excitation mode graviton) 
decaying to a pair of SM Higgses with both the Higgs decaying to a pair of b-quarks in the context 
of the Randall-Sundrum model [87] . The 95% exclusion limits set by the ATLAS collaboration are 
sensitive for heavy resonances (mass > 500 GeV) only, and thus is not directly applicable to the 
parameter space of interest, so we do not consider them in our analysis. 

3.1.4 Search for H/A with final states 

The coupling of the MSSM heavy Higgs {H) and pseudoscalar Higgs (A) with down type fermion 
fd (say bottom quark, tau lepton) is proportional to cos a/cos/3 and tan/3 respectively. So, for 
a fixed value of Higgs mixing angle a, both the couplings Hfdfd and Afdfd increases with tan/3. 
Thus for large values of tan/3 (say > 10), both H and A dominantly decays to bb (~ 90%) and 
r+r“ (~ 10%), resulting into strong suppression in all other decay modes [20]. The production 
of the MSSM heavy Higgses is also primarily controlled by tan/3. The Higges are dominantly 
produced via gluon-gluon fusion process and associated production with b-quarks i.e., 66$ ($= 
H,A). All other production processes like VBF, associated production with gauge bosons, associated 
production with top quarks are dominantly suppressed for the heavy Higgses H/A due to small 
HVV/AVV and Htt/Ati couplings [20]. The ATLAS and CMS collaborations at the LHC have 
studied the signatures of these heavy Higgses [H/A) produced via ggF and b-quark associated 
production processes, and decays to pair of r-leptons [54,55]. No excesses are observed over the 
SM backgrounds, and thus model independent bounds are placed on the production cross section 
times branching ratio a x Br($ —)• for different values of with $ = H/A. In Fig. 9(a) and 

Fig. 9(b), we display our scanned points along with the upper limit on production cross section by 
the ATLAS and CMS collaborations for two dominant production modes, left panel corresponds to 
the ggF process while right panel represents the associated production process 66$, $ = 77, A. In 
order to compare our findings with the ATLAS [54] and CMS [55] data, in Fig. 9(a) and Fig. 9(b) 
we show the distributions of cr x Br($ —)• combining the contributions of both H and A for 

the two different production mechanism. 

It is evident from Fig. 9 that the LHC data from the H —)■ channel has significant impact 

on the 2(7 allowed parameter space. In fact, the entire region with tan /3 > 20 and significant 
amount of parameter space with tan/3 > 10 are excluded when the Higgs is produced via 66$ 
process. However, sensitivity of this channel with ggF production process is negligible. For a better 
understanding of the sensitivity of channel, we display the scanned dataset with explicit tan /3 
dependence, regions with 30 < tan/3 < 50, 20 < tan/3 < 30, 10 < tan/3 < 20 and 1 < tan/3 < 10 
are shown in pink (cross), green (triangle), blue (circle) and violet (square) respectively. One can 
now clearly see that the regions with large tan j3 implies large coupling with 6 and rs and so large 
a X Br, and thus more stringent constraint by the LHC data. 
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Figure 9: Scatter plot in Mh — [cr x Br($ —)■ T+r“)] plane where ^ = H, A produced via (a) gluon 
fusion and (b) in association with b-quarks. The solid red (black dashed) line represents the observed 
upper limits on a x Br($ —)• at 95% C.L. by the ATLAS (CMS) collaboration using LHC-8 

data [ 54 , 55]. Combined contributions of a x Br(t^ —)■ r+r”) of both H and A for the 2a allowed 
parameter space with 30 < tan /3 < 50, 20 < tan /3 < 30, 10 < tan /3 < 20 and 1 < tan /3 < 10 are 
shown in pink (cross), green (triangle), blue (circle) and violet (square) respectively. 


3.1.5 Search for A with Zh final states 

Similar to the case where we discussed the possibility of heavy resonances decaying to a pair of 
MSSM lightest Higgs bosons (see Sec. 3.1.3), one can also consider the possibility of having h 
in association with Z boson from the decay of pseudoscalar Higgs boson A when Ma is greater 
than {Mh + Mz). However, the decay rate r(A —)• Zh) is appreciable only at very low values of 
tan/3 (< 10) and for Ma below the tt threshold i.e., approx. 350 GeV. Once the tt decay opens 
up, it dominates for all values Ma > 350 GeV and tan/3 < 10 [20]. The GMS collaboration 
at the LHC has searched for the heavy pseudoscalar Higgs bosons when it decay to a Z boson 
and a light Higgs boson (h). The final state includes two opposite sign leptons (from Z decay) 
and two b-quarks (from h decay). Thus, one can fully reconstruct the mass of A using the four 
momentum information of the final state leptons and b-jets. In fact three clear distinguishable 
resonance peaks, at Mz, Mh and Ma, are expected for the signal events. Both the ATLAS and 
CMS collaborations have performed the search for the heavy pseudoscalar Higgs bosons at the 8 
TeV run of LHC [60,61]. No significant excesses over SM background are observed in the ATLAS 
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Figure 10: faj Left: Scatter plot in Ma — [ct x Br(A —)• Zh —)■ plane assuming gluon-gluon 

fusion production process for 2a allowed parameter space (represented by magenta points). The 
solid red (blue dashed) line represents the observed upper limits on a x Br(A —)• Zh —)■ at 

95% C.L. by the CMS collaboration using LHC-8 data [61] with narrow width (Ta = 30 GeV) 
assumption, (b) Middle: Scatter plot in Ma — [cj x Br(A —)> Zh) x Br(h —)• bb)] plane for ggF 
production mode. The red solid line indicates the observed upper limits on a x Br(A —)■ Zh) x 
Br(h —7- bb) at the 95% C.L. by the ATLAS collaboration using LHC 8 TeV data [60]. (c) Right: 
Scatter plot in Ma — [ax Br(A —)• Zh) x Br(h —)• plane for ggF production mode. The red 

solid line shows the observed upper limits on a x Br(A —)• Zh) x Br(h —?• at the 95% C.L. by 

the ATLAS collaboration using LHC 8 TeV data [60] . 


and CMS data. Thus model independent 95% C.L. upper limits are imposed on the production 
cross section of A times Br(A —)■ Zh). In Fig. 10(a) and (b), we present the 2cr allowed scanned 
data points along with the CMS and ATLAS exclusion limits respectively. The solid red exclusion 
line in Fig. 10(a) assume narrow-width resonance, while the blue dashed line consider the 30 GeV 
decay width of A boson. The ATLAS collaboration has also performed their study with narrow 
width approximation, however they have analyzed two possible decay modes of the Higgs boson, 
namely bb and t^t~. In Fig. 10(b) we display the allowed parameter space with bb hnal state, 
while Fig. 10(c) shows the same but with final state. As we already mentioned, this decay 

mode plays an important role only in the region with small tan/3 and low Ma- We find that in 
the parameter space of our interest, ATLAS and CMS data are not sensitive enough to impose any 
additional constraints. We expect these heavy pseudoscalar Higgses will be probed with improved 
measurement of these decay modes in the run-H of LHC. 
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3.2 Charged Higgs boson searches 


3.2.1 Search for with tv and cs final states 

The SM particle content does not include a charged scalar particle, however models with additional 
Higgs doublets predict charged Higgs bosons {H^). Thus the discovery of a charged scalar particle 
is a clear signature of the BSM physics. The large electron positron collider (LEP) searched for the 
charged Higgs bosons with center-of-mass energy ^/s = 209 GeV. LEP did not found any signal of 
charged Higgs boson which then leads to put bounds on mass of the charged Higgs boson > 

78.6 GeV [88]. The production and decay of primarily depend on and top quark mass, 

more precisely on whether Mfj± < Mtop or M^± > Mtop. If the charged Higgs boson is lighter 
than the top quark, i.e., Mh± < {mt — mb), then is mostly produced from the tt process. The 
decay of depends on the coupling of the charged Higgs boson with the fermions which is mainly 
controlled by tan/3. The coupling is large for very low and very large values of tan/3 and small 
for intermediate values. Thus, for light enough it is primarily produced from the tt process 
and dominantly decays into TVr final states. It is to be noted that light charged Higgs bosons 
{Mjj± < Mt) can also decay to a charm and a anti-strange quark, and thus in certain regions of 
parameter space one can expect some competition between the cs and tVt- decay modes. Now, if 
the is heavy i.e., > Mtop, then they are mainly produced in associated production with a 

top quark i.e., pp —>■ tH^ + X. For small tan/3, exclusively decays to a top and bottom quark, 
however for large values of tan/3 the decay of to tVj- is not negligible i.e., Br(iL^ —?■ t^Vt) ~ 
10 % [ 20 ]. 

Both GMS and ATLAS collaborations have searched for the charged Higgs bosons at the 8 TeV 
run of LHG using the top-quark pair production process and associated production of with a top 
quark. Search for the light charged Higgs bosons decaying to a r and a tau-neutrino (vr), and/or 
to a charm and strange quark are also presented by the ATLAS and GMS collaborations [62-65]. 
The results seem to agree with the SM predictions and non observation of the excesses leads to 95% 
C.L. exclusion limits on the production cross section times branching ratios for different values of 
Mh±. 

Now, we find that all the points of the scanned data set corresponds to Mf{± > 200 GeV and 
thereby the decay t —>■ bH^ is kinematically forbidden, and so ATLAS and CMS bounds on Mfj± 
using the tt sample have no effect on the parameter space of interest. However, we do impose the 
ATLAS and CMS bounds on Mfj± considering the production of charged Higgs boson in association 
to a top and bottom quark with decaying to t^Vt- In Fig. 11(a), we superimpose the ATLAS 
and CMS bounds on the points corresponding to our data set. We check that the bounds using the 
cs mode are not applicable. From Fig. 11(a) one can infer that one/two orders of improvement in 
cross section measurement might help to probe the parameter space of interest. 
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Figure 11: (a) Left: Scatter plot in Mfj± —[a —)■ T^u-r)] plane for ggF production. Magenta 

points represent 2a allowed regions. The solid red (blue dashed) line represents the observed upper 
limits on a X Br(H^ —)■ t^Vt) at 95% C.L. by the ATLAS [62] (CMS ) collaboration using LHC 8 
TeV data, (b) Right: Scatter plot in Mfj± — [ax Br(H=*= —>■ tb)] plane for ggF production. The solid 
red line shows the observed upper limits on a x Br(H=*= —)■ tb) at 95% C.L. by the CMS collaboration 
using LHC 8 TeV data. 


3.2.2 Search for H^ with th final states 

As we have already mentioned, the decay of H^ to a top and bottom quark dominates in the regions 
with Mij± > Mtop. They are primarily produced via gg —?■ tbH^ process. However processes like 
qq —7- H'^H~, associated production with neutral Higgses can give small contributions to the 
dominant tbH^ process. The decay of H^ is completely controlled by the single parameter tan/3. 
For small values of tan/3 it is the dominant decay mode (branching ratio is close to unity), however 
at high tan/3 the branching ratio decreases slightly Br(Ff^ —>■ tb) ~ 90% while Br(Fr^ —)• ~ 

10%. The CMS collaboration have searched for the charged Higgs bosons with H^ decaying to a top 
and bottom quark [66]. The process under consideration looks like gg —?■ H^tb —?■ {iv^bb){i'vpb)b, 
with £, F being an electron or a muon. The search is performed with 19.7 fb“^ of data at ^/s = 8 
TeV. No evidence for a charged Higgs signal is found and thus upper limits on the production rates 
are placed for H^ masses in the range of 180-600 GeV. We calculate the gg —)• tbH^ production 
cross section for all our valid points, and then impose the CMS bounds. We present our results in 
Fig. 11(b), where the solid red line indicates the CMS exclusion limit. It is evident from the figure 


24 















that tb mode is not sensitive enough to put any strong bound on the parameter space. 


4 SUSY Higgs : Future limits 

4.1 Heavy Higgs search 

In the last section, we discussed the limits on the allowed parameter space of the MSSM Higgs sector 
coming from the direct searches at LHC run-I. We show that significant amount of parameter space 
with additional light Higgses are still allowed by the LHC data. Additionally, we also found that 
the production of heavy Higgs bosons in association with the bottom quarks with HjA decaying to 
T^T~ provides the most stringent bound on our parameter space and the regions with large values 
of tan/3 (> 20) are excluded by the decay mode. In this section, we discuss the sensitivity 

of HL-LHC for the heavy Higgses through the following decay modes: H —)• hh, H ^ tt, A ^ Zh 
and H/A ^ Note that, these channels^ are expected to be most sensitive at the regions 

with tan /3 < 20. Additionally, we also consider the decay H —)• ZZ following the ATLAS and CMS 
analysis [94,95] at 14 TeV run of LHC with 3000 fb~^ of integrated luminosity. 

4.1.1 Search for H with M final states 

Both the ATLAS and CMS collaborations have looked for the heavy Higgs bosons which are pro¬ 
duced via ggF process and decayed into ZZ {Z —)■£/’,£ = e,fj.) at 14 TeV HL-LHC [94,95]. In 
their analyses, the ATLAS collaboration has considered the decay width of H to be the same as 
that of a SM-like Higgs boson for a given mass, while the CMS collaboration calculated the width 
of H assuming tan/3 = 1 and cos(/3 — a) = - 0.06. The expected 95% C.L. exclusion limits on ax 
Br(Lf —7- ZZ —)• At) at £ = 3000 fb~^ by the ATLAS and CMS collaborations are quite similar 
as shown in Fig. 12 by red solid and blue dashed line respectively. We calculate the quantity ax 
Br(Lf —)> ZZ —)• 4£) for all the points in our scanned data set and then compare our results with 
ATLAS and CMS predictions [94,95]. In Fig. 12 we overlay the experimental predictions with our 
scanned data set shown in magenta triangles. We find that for most of the parameter space points, 
even for the points with large production cross sections, the branching ratio Br(H —>■ ZZ) is very 
small (~ 10“^ to 10“®)^, which results into smaller values of ax Br(H —?■ ZZ —)• 4£). We thus find 
that most of the parameter space points are beyond the reach of HL-LHC for H ^ 4£ hnal state. 
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Figure 12: Scatter plot in Mh — [ax Br(H —)• ZZ —)• plane assuming gluon-gluon fusion 

production process. Magenta colored points represents for 2a allowed parameter space. The solid 
red line (blue dashed line) represents the expected upper limits on a x Br(H -^TIL ^ at 

95% C.L. by the ATLAS (CMS) collaboration [94, 95] for I 4 TeV LHC with £ = 3000 fb“^. 

4.1.2 Search for pseudoscalar A with £~^i~bb final states 

In the regions below the tt threshold and small tan/3 the pseudoscalar Higgs A decays to Zh with 
an appreciable amount. An interesting feature of this decay mode with Z —?■ and h ^ bb 

is that one can fully reconstruct the mass A using the four momentum of the leptons and b-jets. 
ATLAS and CMS collaborations have analysed the sensitivity of this channel at the HL-LHC via 
ggF process [94,95]. In order to present the expected 95% C.L. exclusion limits ATLAS has assumed 
a narrow width approximation (i.e., width of A is much smaller than the experimental resolution) 
while CMS has calculated the width of A by assuming tan/3 = 1 and cos(/3 — a) = - 0.06. In 
Fig. 13, we show the distribution of cr x Br(-^ Zh —)• £ibb) for all the points in the scanned data 
set, and then overlay the 95 % C.L. upper limits by the ATLAS (red solid line) and CMS (blue 
dashed line) collaborations [94,95] at 14 TeV LHC with £ = 3000 fb~^. We find that the ATLAS 
limits are more stronger (by almost one order) than the CMS, although the reason behind this is 
apparently not clear. From the Fig. 13 it is clear that only a very small region of the parameter 
space will be excluded by the HL-LHC data, in fact to probe the remaining parameter space few 

■^Reach of heavy Higgs bosons have been discussed in the literature in the context of 8/14 TeV run of LHC, for 
details see Refs. [89-93]. 

® The result is consistent with the alignment limit. 
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Figure 13: Scatter plot in Ma — [ct x Br(A —)■ Zh —>• £+£“bb)] plane assuming gluon-gluon fusion 
production process. Magenta colored points represents for 2a allowed parameter space. The solid 
red line (blue dashed line) represents the expected upper limits on a x Br(A —)• Zb —>■ at 

95% C.L. by the ATLAS (CMS) collaboration for If TeV LHC with C = 3000 fb“^ [Of, 95]. 


orders of magnitude improvement in cross section measurement is required. 

4.1.3 Search for H with di-higgs {H —)• /i/i —)• 6677 ) final states 

In this subsection we discuss the possibility of observing the heavy CP even Higgs boson [H) at 
HL-LHC with its decay to a pair of SM-like Higgses {h) with one Higgs decaying to bh and other to 
77 modes. Single H production cross section can be up to two orders of magnitude larger compared 
to the direct h pair production cross section (see Table 1 of Ref. [96])® depending on the choice of 
model parameters and it can also have non-trivial effects on the self coupling measurement of the 
125 GeV Higgs [96]. 

In the MSSM, the production cross section of H and its decay to a pair of SM-like Higgses 
crucially depends on the SUSY parameter space, mainly Ma and tan [3. More precisely, below the 
tt threshold (350 GeV), the decay rate —)• hh) is substantial only for smaller values of tan/3. 

The most dominant production mechanism of H is the ggF process although, for large or moderate 
tan/3, the bottom quark annihilation to H {bh —)• H) cross section can be substantial [97]. The 
final state signature, depending on the decay of the SM-like Higgses, includes, for example, bbbb, 

®The SM Higgs pair production cross section at next to leading order is about 34 fb at 14 TeV LHC for Mn = 125 
GeV [98]. 
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bbT~^T~, bbW^W~, 6677 etc. Among all these possibilities, bbbb final state has the largest cross 
section. However, due to enormous QCD background this is one of the most challenging scenarios 
to be observed at the LHC^. On the other hand, even though the branching ratio for 6677 channel 
is very small (about 0.27 %), it is the most promising channel due to large photon identification 
efficiency and very good resolution in the photon energy measurement. Hence using the di-photon 
invariant mass distribution one can easily reconstruct the Higgs boson mass and at the same time 
separate the signal from the SM background®. 
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Figure 14: Sensitivity of di-Higgs final state in Ma — tanfi plane at HL-LHC from bb'j'j channel 
with L =3000 fb“^. The magenta coloured points are 2a allowed points from global analysis and 
the black circled points are expected to be probed at HL-LHC (see Sec. 4-1.3 for details). 


A detailed signal-background analysis of heavy Higgs production and its decay to a pair of 125 
GeV Higgs in 6677 channel has been already performed in Ref. [96]. Here we use their results to 
constrain our parameter space®. Events with two 6 -jets, two photons and no isolated leptons are 
selected after imposing the basic selection cuts following the ATLAS collaboration [101]. Using the 
four momentum information of the bb and 77 system, the invariant mass of the heavy Higgs boson 
can be reconstructed assuming is Mh ± 50 GeV. 

^ It has been shown that the jet substructure technique can be very useful to separate the signal events from the 
large QCD and electroweak backgrounds [99]. 

®In this case, the dominant backgrounds are tih and the direct Higgs pair production (hh). Note that, the 
production cross section for both the Higgs pair and tih processes depend on the MSSM parameters, however, here 
we assume SM cross sections for hh and tih processes. 

®The analysis for 6677 channel was hrst introduced in the Ref. [100] 
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Let us now estimate the sensitivity of the HL-LHC to probe the parameter space in the 6677 
channel. In Ref. [96], authors considered few representative benchmark points with heavy Higgs 
mass Mh in the range of 275 - 600 GeV and performed a detailed collider analysis. We should 
note, as also pointed out in Ref. [96], the cut efficiencies do not change radically in the entire region 
of 275 - 600 GeV. The number of background events (Nb) and the cut efficiencies as a function 
of M[{ are taken from Ref. [96]. We estimate the number of signal events {Ns) for all the points 
in the scanned data set by multiplying the cut efficiencies with the production cross sections, and 
then calculate the signal significance S = Ns/y/Ns- In Fig. 14, we present the sensitivity of the 
bb'y'y channel at HL-LHC where magenta triangles corresponds to our scanned data set, while the 
black circles represent the points with S > 2a. In other words, these are the points (black circle) 
with low tan/3 (< 10) that are expected to be probed at the HL-LHC. The lack of sensitivity of 
HL-LHC in the regions of parameter space with tan /3 > 10 is due to the fact that both the ggF 
production cross section and Br(77 —?■ hh) decreases with the increase of tan /3 forcing this portion 
of parameter space to go beyond the reach of HL-LHC. From the Fig. 14, one can note that future 
run of LHC might not be able to exclude the entire region below Mh ^ 425 GeV and tan /3 ^ 8 . 

4.1.4 Search for H/A with ti final state 

Above the kinematic threshold of tt (~ 350 GeV), the decay of H and A into tt pair opens up. In 
fact, for Ma > 350 GeV and low to moderate values of tan/3, it is indeed the dominant decay mode 
of H/A. Besides, due to large top-Higgs Yukawa coupling, the production cross section of the heavy 
Higgs H via ggF process also becomes large in the low tan (3 regime. Hence, assuming the narrow 
width approximation, one can expect to observe a resonance peak at Mh/a the tt invariant mass 
distribution. However, the main drawback of such a bump hunting procedure is that it is extremely 
difficult to extract the ti resonance peak from the huge SM ti continuum background. Here, we 
perform a detailed signal-background analysis in the context of HL-LHC and study the sensitivity 
of HL-LHC to probe the region of parameter space of our interest. 

We analyze the production of heavy Higgses H/A via ggF process and their decay to 
where one top decays leptonically and one decays hadronically. So, the final state includes one 
isolated lepton (electron or muon), at least four jets among them at least two are 6 -jets and missing 
transverse energy (^-p). We use PYTHIA (version 6.4.28) [ 86 ] to generate both the signal and 
dominant SM background, e.g., ti events. Electrons are selected with pp > 20 GeV and |ry| < 2.47, 
while we choose muons with pT > 20 GeV and |r/| < 2.4. We select jets with pT > 30 GeV and 
\r]\ < 2.8. A jet is called a “b-jet” if the angular separation AR = sj (Ar/)^ -|- {A(j))^ in the r] — c/ 
plane between the jet and a B-hadron is less than 0.2. Following the ATLAS collaboration, we 
also assume a flat 70% b-tagging efficiency for the b-jets [102]. In our analysis, we also implement 

H and A are almost degenerate, one might not be able to distinguish them in tt invariant mass distribution. 
For this reason, here we consider the production of the heavy Higgses H/A simultaneously via the ggF process. 
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Figure 15: Distribution of tt invariant mass for some representative SUSY benchmark points 
with Mh= 400 GeV (green), 500 GeV (blue), 600 GeV (red) and SM tt background (black). The 
distribution has been drawn before applying the cut on invariant mass (Ch) in signal region SR- 
loose (see text for details) for integrated luminosity C = 3000 fb“^. We consider NNLO+NNLL 
cross section for SM tt production at If TeV LHG (966 pb) and for signal the no of events has 
been estimated assuming single H production cross section to be 1 pb with Br{H tt) = 100%. 


the methodology of lepton isolation and lepton-jet identification following the ATLAS study [103]. 
Although it is possible to reconstruct leptonically decaying top quark within a quadratic ambiguity, 
for simplicity we use the parton level four momentum information of the neutrino to reconstruct 
the H mass. 

We are now in a position to describe the details of our simulation procedure as well as the kine¬ 
matic selection cuts for our signal and the backgrounds. We choose three representative benchmark 
points with Mh = 400, 500 and 600 GeV, and define three signal regions (SR), namely SR-loose, 
SR-medium, SR-tight, depending on the choice of our selection cuts. The SR-loose is defined based 
on the following set of cuts: 

• Cl: Events must contain at least 4 jets with px > 50 GeV. 

• C2: Among the four selected jets, two are b — tagged. 

• C3: Events containing one isolated lepton with pT > 30 GeV are selected. 

• C4: Missing energy (^^rp) > 50 GeV. 

• C5: Select events with tt invariant mass between Mh^ 25 GeV. 
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We denote the signal region SR-medium with the same set of cuts as that of SR-loose but with 
> 100 GeV. In order to probe higher values of Mh, we define the signal region SR-tight where 
the cuts on pT of the jets and are stronger. For example, for SR-tight, events are selected with 
Pt of the first two leading jets greater than 100 GeV while the missing transverse energy ^rp > 100 
GeV. 


Ghannel 

Number of Events at 3000 fb ^ 

Mh = 400 GeV 

Mh = 500 GeV 

Mh = 600 GeV 

Signal 

tt 

Signal 

tt 

Signal 

tt 

SR — loose 

1268 

104612 

9658 

420572 

26842 

563452 

SR — medium 

8 

1741 

1584 

69232 

9656 

194698 

SR — tight 

- 

- 

4 

637 

2296 

44894 


Table 7; Number of signal (H production via ggF) and SM tt events after applying the cuts in SR-Loose, 
SR-Medium and SR-Tight signal region at 14 TeV LHC with C = 3000 fb~^. For tt, we use NNLO cross 
section [104]- For the signal estimation we present the number assuming <7{pp —>■ H)nlo x Br{H tt) = 
1 pb. 


In Fig. 15, we display the tt invariant mass distribution for the three representative benchmark 
points with Mh= 400 GeV (green), 500 GeV (blue), 600 GeV (red), and also overlay the same for 
the SM ti background (black line). We show the invariant mass distribution for the signal region 
SR-loose with an integrated £ = 3000 fb“^. To estimate the number of background events, we use 
the ti cross section at next-to next-to leading order (NNLO) ~ [1^*4]. For simplicity, 

we assume that for all the benchmark points Br{H —>■ ti) is 100% and the next-to leading order 
production cross section o'{pp —>■ i7) is 1 pb. Note that, we make such a conservative choice just to 
keep our analysis simple, one can easily scale our numbers with the actual production cross sections 
and branching ratios. Although one can see clear resonance peaks at the Mh masses (see Fig. 15), 
enormous SM ti background makes it very challenging to observe a clear signal of heavy Higgs 
in the ti invariant mass distribution over the SM background. To be more precise, we count the 
number of signal and background events for three signal regions SR-loose, SR-medium and SR-tight 
after G5, and then calculate the statistical significance S = Ns/y/Ns- In Table 7, we present the 
number of signal and background events for three benchmark points and the SM ti backgrounds at 
the 14 TeV run of LHG with £ = 3000 fb“^. From Table 7 one can see that Ns/Nb ratio is very 
small for all the benchmark points irrespective of the signal regions. We find that the statistical 
significances S for Mh = 600 GeV are 36, 22 and 11 for the three signal regions SR-loose, SR- 
medium and SR-tight respectively. However, even with 5% systematic uncertainty these numbers 
reduces to 0.95, 0.99, 1.02 respectively. 

In Fig. 16 we show the distribution of cr x Br[H, A ^ ti) assuming production of both H and 
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Figure 16: Scatter plot in Mh — [cr x Br(H, A —)• tt)] plane assuming both H and A production via 
ggF at 14 Te V LHC. Magenta colored points represents for 2a allowed parameter space. Red line 
corresponds to a x Br{H, tt) = 1 pb. 


A via ggF for all the points in our scanned data set. We find that a x Br{H, A ^ tt) lies mostly in 
the region 0.5 - 0.0001 pb. The red solid line in Fig. 16 represents a x Br(H, A ^ tt) = 1 pb. This 
red solid line indicates that for most the points in our scanned data set the quantity cross section 
times Br(H, A —t- ft) is 10 - 1000 times smaller, and thus the the numbers in Table 7 represent 
a too much optimistic scenario. In other words, even HL-LHC might not be sensitive enough to 
probe such a region of parameter space. 

Recently, it has been shown that using angular cuts, the signal significance can be improved 
[105]. However, one should note that the inclusion of systematic uncertainties may change the 
significance drastically and it is not possible to predict the reach of the ti channel at the HL-LHC 
without the precise knowledge of the systematic uncertainties. We end this section by mentioning 
that the observation of a heavy Higgs in H ^ tt channel at the HL-LHC is really a challenging 
task and it needs special attention and more detailed studies. 

4.1.5 Search for H/A with final states 

In Sec. 3.1.4 we study the impact of the H/A —)• channel on the parameter space of interest 

using the LHC-8 data and find that regions with large tan/3 (>20) and Ma up to 600 GeV are 
already excluded. In this section, we discuss the sensitivity of the 14 TeV run of LHC to probe the 
remaining allowed region of parameter space via the channel. In order to perform a detailed 
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signal-background analysis, we follow the ATLAS simulation with LHC-8 data [54]. We assume 
that H/A are produced via the b-quark associated production process and decays to r+T“. The 
ATLAS analysis [54] have shown that for Ma between 300 and 600 GeV, the final state with one 
leptonically decaying r (riep) and one hadronically decaying r (rhad) gives the best sensitivity (see 
Fig. 9(b) of Ref. [54]), and so in our analysis we consider the H/A ^ ^lepThad signature only. We 
identify taus through their hadronic decays and demand that the candidate jet must lie within 
\rj\ < 2.5 with pT > 30 GeV. The jet must contain one or three charged tracks with [r/trackl < 2.5 
and highest track pT > 3 GeV. Moreover, in order to ensure proper charge track isolation, we also 
require that there are no other charged tracks with pT > 1 GeV inside the candidate jet. 

We select events with at least one lepton (electron/muon) with pT > 50 GeV and an oppositely 
charged r-hadron with pT > 50 GeV. We further demand that events do not contain any additional 
electrons or muons, and the transverse momentum difference Ap^ = (rhad)~P t( lepton) is greater 
than 50 GeV. Additionally, following Ref. [54], we also demand the sum of the azimuthal angles 
Y//Ac/ = A</)(rhad) ^t) + ^'/’(Repj^x) than 3.3, and the hadronic and leptonic r-decays 

satisfy A(/)(riep, Thad) > 2.4. Besides, events with at least one b-tagged jet with pT > 50 GeV are 
also selected. Finally, we demand that the reconstructed di-tau invariant mass (rurr) is within 
Af$± 30 GeV (‘k = H/A). The dominant SM background processes that can mimic the signal 
are Z -l-jets, W -l-jets, QGD multi-jets and tt. We analyze the most dominant SM backgrounds, 
namely Z -|-jets and tt, while we check that QGD multi-jet contribution is negligibly small. In 
order to accommodate other subdominant backgrounds, we consider 50% enhancement of the total 
estimated background. 

The signal events are generated using PYTHIA (version 6.4.28) [86] assuming H/A are produced 
via the b-quark associated production process. Here we concentrate on the region 200 < Ma < 600 
GeV and split the entire Ma region in steps of 25 GeV. For each event we then calculate the 
di-tau invariant mass using the “collinear approximation technique” with the assumption that the 
r-lepton and all its decay products are collinear [106] and in an event is entirely due to the 
neutrinos^^. After performing a bin-wise signal/background analysis, we thus obtain a 95% G.L. 
expected exclusion limit on the production cross sections for the signal events assuming BR(H /A —)■ 
r+r“) = 100%. We calculate the production cross sections for all the points corresponding to our 
scanned data set, and then impose the estimated exclusion limits on the production cross sections 
for different values of Ma (or Mh)- In Fig. 17(a), we display the distribution of the quantity 
a X Br{H/A —?■ Mt~) when H/A are produced via the b-associated process. Similar to Fig. 9, we 
here also display the scanned data set with different tan/3 dependences. For example, the regions 
with 20 < tan/3 < 50, 10 < tan/3 < 20, 5 < tan^d < 10 and 1 < tan/3 < 5 are shown in red (cross), 
pink (triangle), blue (circle) and green (square) respectively. The solid red and black lines indicate 

Another technique often used by the experimental collaborations are the Missing Mass Calculator (MMC) method 
which allows a better reconstruction of the tt invariant mass distribution [106]. In this paper, however, we restrict 
ourselves to the “collinear approximation technique”. 
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Figure 17: Left: Scatter plot in Mh — [ct x Br(<J> — )• r+r”)] plane where ^ = H, A produced in 
association with b-quarks. The entire tan /3 region is splitted into different pieces, the regions with 
20 < tan/3 < 50, 10 < tan/3 < 20, 5 < tan/3 < 10 and 1 < tan/3 < 5 are shown in red (cross), pink 
(triangle), blue (circle) and green (square) respectively. The solid red (black) line represents the 
expected upper limits on a x Br($ —)■ Mt~) at the HL-LHC with C =3000 fb“^. Right: Sensitivity 
of H/A —)• T'^T~ final states in Ma — tanfd plane. The magenta (plus) points are the 2a allowed 
points corresponding to our scanned data set, while the blue (square) and green (cross) points are 
expected to be probed at HL-LHC considering 20% and 10% systematic uncertainties respectively. 


the expected sensitivity in the HL-LHC with 10% and 20% systematic uncertainties respectively. 

We have already mentioned that in the MSSM, the production cross sections and branching 
ratios of H and/or A crucially depend on two important parameters Ma and tan/3. Thus the 
expected sensitivity of the 14 TeV LHC on the H/A production cross sections can be translated in 
the Ma — tan/3 plane. In Fig. 17(b) we display such possibility where the magenta (plus) points 
denote our scanned data set while the blue (squared) and green (crossed) points represent the 
regions of parameter space in the Ma — tan /3 plane which can be probed at the high luminosity 
run of LHC with 3000 fb“^ of data assuming 20% and 10% systematic uncertainties respectively. 

Thus, combining the existing 8 TeV direct search bounds (see Fig. 9) and the estimated HL- 
LHC bounds, one can infer that the regions with tan/3 > 20 are already excluded by LHC-8 data 
while the regions with tan/3 down to 8 with any values of Ma can be probed at the HL-LHC. 
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5 Conclusion 


A scalar particle with mass close to 125 GeV has been discovered at the LHC. Measurements of 
spin-parity and various couplings seem to agree with the SM expectations. In the minimal extension 
of SM, namely MSSM, one can identify the observed 125 GeV Higgs boson as the lightest Higgs 
boson among the five MSSM Higgses h, H, A and H^. In the MSSM, the couplings of h with 
the SM gauge bosons {WjZ) are proportional to sin(/l — a) where tan/3 is the ratio of vevs of 
two Higgs doublets while a is the Higgs mixing angle. Precise measurements of various couplings 
of the observed Higgs boson with WjZ bosons by the ATLAS and CMS collaborations imply 
sin(/3 — a) ~ 1. At the tree level, the Higgs mixing angle can be derived using pseudoscalar Higgs 
mass parameter Ma and tan/3, however if we include radiative corrections, a becomes a non-trivial 
function of various SUSY parameters. The ATLAS and CMS collaborations have also searched 
for the heavy Higgs bosons (iL, A, H^), however absence of any signal puts strong bounds on the 
masses/branchings of these heavy Higgses. One can easily satisfy the current LHC data by taking 
a —7- 0, and /3 —)• 7r/2 with Ma >> which is generally known as the decoupling limit of MSSM. 
In this limit, the lightest MSSM Higgs boson behaves exactly like SM Higgs and masses of other 
Higgs bosons (iL, A, are pushed well above LHC reach. In this paper we study the possibility 
of having light additional MSSM Higgs bosons, preferably below 600 GeV, with moderate Higgs 
mixing angle a, being consistent with the SM Higgs data and also direct search limits on the MSSM 
heavy Higgses. 

We restrict ourselves to the 19 dimensional pMSSM framework and scan the parameters that 
are relevant for the MSSM Higgs sector. We perform a global fit analysis using most updated data 
(till December 2014) from the LHC and Tevatron experiments and also consider the flavor physics 
constraints. The region with Ma < 350 GeV and tan/3 > 25 are excluded by the Br(Bg —)■ 
while Ma < 350 GeV with tan/3 < 8 is not favoured by the Br(b —)■ sy) constraint. The regions 
with large Ma value (> 400 GeV) are not much constrained by the data. An interesting point to 
note that regions with 200 < Ma < 400 can have moderate values of the Higgs mixing angle a = 
0.2, while for relatively large values of Ma, ot can be as large as ~ 0.8 with small tan/3. Thus, one 
is not always forced to be in the decoupling limit to comply with the LHC data and light additional 
Higgses (Ma) are still allowed by the current data. Moreover, 10 - 20% deviations from the SM 
expectations are also observed for various Higgs signal strength variables. 

We next study the impact of current bounds on the MSSM heavy Higgs boson {H, A and H^ ) 
masses and couplings from the direct search at the LHC. We analyze the following decay modes of 
the MSSM heavy Higgses; H —^ 77 , H —^ WW, H ^ hh ^ bbbb, H ^ hh ^ btr/^, H/A^ , 
A Zh, and —)■ tb. As we have already mentioned, most of the regions of parameter 

space satisfy the alignment limit (/3 — a ~ 7r/2), the 7-1-8 TeV LHC data on H WW channel is 
not sensitive enough to constrain the parameter space of our interest. Except H/A —)■ , LHC 

bounds on the production cross section times branching ratios for other decay modes of H, A and 
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are also one/two orders of magnitude larger compared to the MSSM expectations. The 7+8 
TeV LHC data on MSSM heavy Higgs boson decay $(= H/A) —)■ put the most stringent 

bound on our parameter space. The entire region with tan /3 > 20 is excluded when the heavy 
Higgses = H/A are produced via process. The reason being large tan/3 (typically > 20) 
implies large Bt:{H/A —)■ r^r“), typically ~ 10%, and large cross section for associated production 
of H/A with bottom quarks, and thus stronger constraint on the quantity a x Hr. 

We also study the prospect of probing the allowed parameter space, mostly low tan /3 region, at 
the high luminosity run of LHC (HL-LHC). We impose the ATLAS and CMS preliminary results 
on if ZZ and A ^ Zh at the HL-LHC. The impact of the future limits from the H ^ ZZ 
and A —)■ Zh channels on the allowed parameter space are found to be marginal at £ = 3000 fb~^. 
Searches via H ^ hh at the HL-LHC can only probe a small part of the parameter space with low 
tan/3. Among several other possible decay modes, H ^ ti dominates in the low tan/3 region for 
Mh > 350 GeV. We perform a dedicated signal-background analysis on the H ^ tt channel by 
choosing few representative MSSM benchmark points, and find that the signal to background ratio 
is very small. A more detailed analysis is required, for example one can use the jet substructure 
technique, spin correlation technique to achieve better sensitivity in this channel. We find that 
so far the combined 7+8 TeV LHC data on the H/A ^ channel has the best sensitivity to 

place constraints on the MSSM parameter space of interest. Thus, one may expect to find stronger 
constraints on the MSSM parameter space from the mode at the HL-LHC. We study this 

possibility and find that combining the existing 8 TeV direct search bounds and the estimated 
HL-LHC bounds, one can infer that the regions with tan/3 > 20 are already excluded by LHC-8 
data while the regions with tan/3 down to 8 with low to moderate values of Ma can be probed at 
the HL-LHC. Below tan /3 < 8, searches via if —)■ /i/i at the HL-LHC can be very important to 
probe the parameter space. 

In summary, we find that regions with low to moderate values of tan/3 with light additional 
Higgses (mass < 600 GeV) remain unconstrained by the current data. Even the high luminosity 
run of LHC may not have enough sensitivity to probe the entire low tan/3 region of parameter 
space. However, the proposed e'’'e“ international linear collider (ILC) will be an ideal machine to 
study this scenario. With the expected accuracies in the determination of various partial decay 
widths of the Higgs boson and also the possibility of producing some of the heavy MSSM Higgses 
directly, one might be able the probe the remaining region of the allowed parameter space with 
^/s = 1000 GeV with higher luminosities at the ILC. 
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